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ABSTRACT 


The infrared vibrational spectra of HOD and DOD dissolved in carbon 
tetrachloride and carbon disulfide were obtained and compared with the 
Same spectra of HOH. The spectra of these three water species were com- 
pared to those of the vapors and the similarities noted. The bending 
mode, 75; of DOD in carbon disulfide was observed. Temperature effects 
on the spectra were also studied. To analyze the data a simple theoret- 
ical vibrational treatment was used providing a computation readily 


adaptable to digital computers. 
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1. Introduction 
Infrared Spectra of HOH 

The infrared spectrum of water has been intensively studied in 
vapor phase, liquid phase, and in dilute solutions of non-polar solvents. 
The vapor phase spectrum has been analyzed by many investigators and in 
considerable detail first by Darling and Dennison. |1| 

The water molecule is well known to have the form of an isosceles 
triangle with equivalent OH bonds, an example of Coy Symmetry. Such a 
molecule possesses three normal modes of vibration, all of which should 
be active in both the infrared and Raman spectra. The H-O-H bond angle 
is reported to be 104.5 + 0.10° with an equilibrium bond distance of 
0.9580A. L-3 The bond angle of DOD is known to be identical. B 

Benedict, Gailar and Plyler present the most recent definitive 
infrared spectral data for HOH, giving the values of the fundamental 
water vapor vibrational frequencies as: Y i= 3656.65 om7!, Y 9 = 1594.56 
en7!, and Diy = 3755619 on™}, | 4) These vibrational modes may be de- 


picted with reasonable accuracy as shown in Figure A. 


7 Vp, % 


Figure A 
Herzberg summarizes previous values of infrared and Raman data for 
the various water molecules HOH, HOD, and DOD. 
Liquid water exhibits the characteristic very broad, intense infra- 


red absorption found for hydroxy containing compounds such as alcohols 


and carboxylic acids. This band is attributed to associated or hydrogen 
bonded species. |6| 
Infrared Spectra of Water in Non-Polar Solvents 

Upon dilution of the OH containing compounds with non-polar solvents 
such as carbon tetrachloride or carbon disulfide, a new, narrow absorp- 
tion band appears at a higher frequency than that of the original associa- 
ted species. 7- | This new absorption is assigned to the Stretching 
frequency of the unassocilated OH bond. 

Infrared spectra of dilute solutions of HOH in carbon tetrachloride 
were studied by Borst et al., who observed two absorptions in the OH 
stretching region at 3702 em7! and 3611 em7!. ho These investigators 
suggested that the 3702 en7! absorption be assigned to 13 and the 3611 
em! band to Y 1- Later Fox and Martin confirmed these observations at 
frequencies of 3705 em~+ and 3614 em! and stated that these absorptions 
were, respectively, the V3 and 7, vibrations of the monomeric HOH 
species in carbon tetrachloride solution. {11 Reisinger has confirmed | 
these observations assigning 3705 em! to V3 and 3613 em! to V1- fa 
None of these investigators reported an observation of the ‘YY. or bend- 
ing frequency. This is to be expected, since in the path lengths used 
for these investigations carbon tetrachloride and carbon disulfide 


absorb intensely in the region near 1595 em7! 


» the vapor phase “Y, fre- 
quency. 

To avoid these solvent absorptions the present investigation utilized 
DOD and HOD whose gas phase V5 frequencies are 1178 em~! and 1402 em7! 
respectively. This yielded an interference free spectrum only for the 


Yo mode of DOD in carbon disulfide, the Y> of HOD being obscured in 


all solvents used. 


Fox and Martin were able to separate graphically a third band in 
the HOH spectrum centered near 3810 em7}, | Noting that this absorp- 
tion appeared to follow the general shape of rotation lines in this 
region they proposed that this absorption was the envelope of ae rota- 
tion lines. Reisinger also noted this third band placing it at 3804 
en7!. 22 Additionally, he separated a fourth absorption in this region 
centered at 3540 cm7!. These third and fourth absorptions were also 
noted by Reisinger in an investigation of a dilute solution of HOH in 
carbon disulfide. [12] The nature of the absorptions appeared identical 
to those of HOH in carbon tetrachloride except for a slight shift of 


1 and 3547 cm7!. Earlier investigators do 


the frequencies to 3796 cm 
not indicate that their studies covered the low frequency region (near 
3500 em7!) so that this fourth band (3540 cm!) was unobserved by them. [11] 

The Yo frequency for HOH in carbon tetrachloride and carbon disul- 
fide would be expected at about 1595 em! (the vapor phase io) since, 
as 1s reported later, the Y9 frequency for DOD in carbon disulfide is 
very little shifted from its vapor phase value of 1178 oe However , 
V2 has not been observed for HOH in carbon tetrachloride or carbon di- 
sulfide. 

Temperature Effects 

Early investigators of the spectra of OH bearing species noted 
that the relative intensities of the OH stretching absorptions were 
functions of temperature. They attributed the effect to temperature 
induced changes in monomer-polymer equilibrium or to solvent effects, 
i.e., a change in the molar absorption coefficient. [8,13 


More recently, however, Slowinski and Claver have noted that for 


pure liquids the intensities of OH stretching absorptions as well as 


C-H, C-Cl and other bond stretching absorptions depend markedly on temper- 
ature. hal These workers show clearly that both the maximum optical 
densities and the integrated intensities of these stretching absorptions 
decrease with rising temperature. Their work was rapidly followed by 
the efforts of other investigators who have extended and refined these 
observations, studying hydrogen bonded equilibria in dilute solutions. 
[15,16 Hughes, Martin, and Coggeshall demonstrate that the intensity 
of the absorbance of the unassociated OH stretching bond decreases lin- 
early with increasing temperature. hs 
Calculations 

The theoretical bases for correlating vibrational frequencies of 
molecules with their parameters (masses and geometries) has been known 
for many years. Wilson, Decius, and Cross have thoroughly delineated 
the methods for these calculations. h7) The matrix methods presented 
by these authors are readily adapted to digital computer technology 
and presented an opportunity to determine if there was a simple theoret- 
ical vibrational treatment consistent with the results of this investi- 
gation. In one simple treatment of a non-linear triatomic molecule 
such aS water, the three fundamental vibrational frequencies are deter- 
mined by four force constants and the molecular parameters. These three 
frequencies are, of course, insufficient information for solution of 
the equations determining the four force constants. However, applica- 
tion of the Born-Oppenheimer approximation permits the assumption that 
the potential function is unchanged by isotopic substitution. This 
allows the six fundamental vibrational frequencies of HOD and DOD to be 
used in determining these same force constants. In the two calculations 


made here the first used the five observed frequencies of HOH and DOD. 


The second used, in addition, the two observed frequencies of HOD for 


a total of seven. 


Ti 


2. Experimental 

Infrared Absorption of the Water Molecules HOH, DOD, and HOD 

The Perkin-Elmer Model 221 recording spectrometer with a sodium 
chloride prism-grating combination was used in the double beam mode of 
operation. The spectrometer was calibrated in the regions of interest 


using the vapor phase spectra of HCl, HBr, NH3, CO,, and HOH as indicated 


9? 
in Table I. The pressures were chosen for the purpose of direct compar- 
ison with the reference values given by the Commission on Molecular 
Spectroscopy. is All gases used except the HOH were Matheson Co. bot- 
tled gases. The HOH spectrum was obtained from atmospheric vapor using 
the spectrometer in single beam operation. The standard deviations from 
the mean of the observed frequencies are given in Table I for each range. 
A mean correction was determined for each range and applied to the ob- 
servations to yield the reported values. 

Connecticut Instrument Corporation 50 mm path length glass cells 
were used with BaF, windows. These cells have a capacity of approxi- 
mately 30 ml. Only glass, Teflon, and BaF, were in contact with the 
Sample in these cells. 

The carbon tetrachloride and carbon disulfide used as solvents 
were Fisher Scientific Co. spectral grade reagents. No additional puri- 
fication was necessary. It was necessary to dry the solvents as com- 
pletely as possible using the disappearance of the 3 and VN; HOH fre- 
quencies as a minimum criterion of dryness. Previous investigators 
have used a variety of techniques for drying, most of which proved un- 
satisfactory when used in this laboratory. 9,12, 19] These included 


drying over CaCl., Na S03, Na,SOy » bubbling with dry air, and distilla- 


Z 


tion. P,05 was a successful drying agent, but its use involves consid- 


12 


TABLE I 


Calabpration Daka 


Standard Frequency (em71) (em!) Pressure (mm) 
HCl 2680-2900 Oig5 120 
HBr 2600-2720 0.5 140 
NH 3 1160-1200 Oy. 3 50 
H90 3745 0.9 * 
CO. 3620 0.9 400 


* Atmospheric vapor 

erable tedious handling. Linde Co. type 4A molecular sieve in bead 
form was used in this investigation, the beads being placed directly 

in a glass stoppered bottle of the solvent. In order to prevent leak- 
age of moisture into the sample, the bottle was placed in a dessicator 
containing P.0s5. No change was observed in the dried solvents or the 
molecular sieve. The solvent spectra were also unchanged over the re- 
gion 4000 to 750 cm! except, of course, for the HOH bands. All trans- 
fers were made uSing a syringe to minimize contact with atmospheric 
moisture. 

Dry carbon tetrachloride with a small amount of added DOD produced 
an observable HOD absorption in as short a time as fifteen minutes when 
placed in the glass-stoppered BaF, cell. However, dry carbon tetra- 
chloride alone did not produce an HOH absorption for a minimum of several 
times as long. It was proposed that exchange of D for H with the glass 
cells may have been responsible for this phenomenon. To prevent con- 
tamination by absorption of atmospheric moisture a few granules of 


molecular sieve were placed in the dry reference cell sample. This 


is 


effectively prevented contamination during a span of several hours, which 
was the longest period observed. All spectra were recorded with an 
identical cell filled with dry solvent inserted in the reference beam 
of the Spectrometer. 

Saturated solutions of HOH and DOD in both solvents were prepared 
by adding excess HOH or DOD to the dried solvent and allowing a minimum 
of 24 hours for equilibration. The DOD was obtained from Bio-Rad Lab- 
oratories and stated to be 99.83 mole percent DOD. 

The various sample solutions were prepared as a given percent of 
saturation by dilution of the saturated solution with dry solvent. 


The equilibrium for the reaction 


2 
: HOD) 
HOH + DOD =? 2 HOD with K EE ||:.0)))) Semen 
[HOH] (DOD) 
is known to be very rapid. The«auilibrium constant has been determined 


as K = 3.96 at “25°C. 20] For a given sample we have 


[yop] + [Hou] + [pop] +c 

where C iS a constant. Noting that 
[Hop]? = K [Hox] [pon] 

it can be seen that the condition for HOD to be a maximum is 
[HoH] = [DoD]. 


So then 


VK. 
[HOD] ax . 2+ VK 


and for K = 3.96 the maximum [HOD]Jis about 50 mole percent, which is 

a ratio of 1:2:1 for the various water species. This maximum HOD con- 
centration exists for an equimolal mixture of HOH and DOD. The values 
in Table II give the observed frequencies for the HOD spectrum in car- 


bon tetrachloride. The notation of Benedict et al. has been adopted 


Frequencies of vibrational absorptions 


TABLE II 


of water (in em!) 


Values given for vapor phase spectra are those of Bene- 


dict et. al. (rounded to the nearest em™1). [4] 


HOH 


Boag 


Tes 


3756 


HOH 


3614 


3705 


HOH 


3603 


3695 


Vapor Phase 
DOD 
2671 
1178 


2788 


Carbon Tetrachloride Solutions 
DOD 


2642 


Zon. 


Carbon Disulfide Solution 
DOD 
2635 
1180 


2744 


15 


of various species 


HOD 


2127 


3707 


HOD 


2692 


3660 


in describing the HOD vibrational bands, namely Vv; = OD stretching, 
V5 = bending, and “Y, = OH stretching vibrations. g 


There is considerable overlap of the V; HOD frequency, 2692 om), 


with those of the DOD spectrum and of the 3 HOD frequency, 3660 en7!, 
with the HOH spectrum. Consequently, the optimum concentration of HOD 
is one that will minimize the interfering bands while allowing sufficient 
intensity of the HOD absorption. A 9:1 misture of HOH and DOD will 
yield approximately an 81:18:1 ratio of HOH:HOD:DOD. This 18:1 ratio 
of HOD to an interfering species is an order of magnitude improvement 
over the previous ratio of 2:1 at maximum mole percent HOD. 
Temperature Effects 

In order to observe the effect of temperature on the spectra 
special cells were constructed. These were 50 mm path length glass 
cells with a coolant jacket as an integral part of the cell. Water 
from a constant temperature bath was used as the regulating fluid. 
Silver chloride windows were held in place onthe ends of the cell body 
by aluminum discs fastened to each other with screws. Internal cell 
temperatures were measured with a thermocouple. Carbon tetrachloride 
solutions were studied over a range from 3°C to 60°C. Due to the high 


vapor pressure encountered in this temperature range, carbon disulfide 


solutions were not studied. 
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3. Results and Discussion 

The results of earlier workers were confirmed upon analysis of the 
spectrum of a dilute solution of HOH in carbon tetrachloride. fro-12] 
This is shown in Figure 1. Two of the fundamental vibrational frequen- 


cles were identified as VY, at 3614 em7! 


and V3 at 3705 em7! in close 
agreement with previously reported values. 11, 12| Easily identified 
also were the broad sidebands centered at 3804 cem-! and 3540 cm™!. These 
sidebands were graphically separated from the spectrum. They were also 
observed by this technique by Reisinger in a solution of HOH in carbon 
disulfide with centers at 3796 cm7! and 3547 cm7! respectively. 9 Fox 
and Martin observed only the 3804 en~! sideband considering it to be the 
envelope of rotation lines of the water molecule. fy] Several possible 
sources for the 3540 cm! band were considered by Reisinger including: 
(1) the possibility of Fermi resonance between 27, and Yy> (2) that 1@ 
is 2V 55 and (3) that it is the envelope of rotation Lines. [12 

Fermi resonance, a quantum mechanical phenomenon, is the result of 
interaction of closely spaced energy levels. In this process these 
levels are spread apart and associated transitions tend to become more 
nearly equal in intensity. A necessary condition for Fermi resonance 
is that the levels involved be of the same symmetry species or share 
at least one species in common. Since 2VY%, and % share the fully 
symmetric species Ay> Fermi resonance is possible if the energy levels 
are sufficiently near one another. 

A determination of the absorption of DOD in carbon tetrachloride 
at varying percentages of saturation at 280C confirmed that the bands 
at 2751 em7! and 2642 em~! should be assigned to the V3 and V1 fre- 


quencies of the DOD molecule. These results are shown in Figure 2. It 


TR 


is observed that the structure of the DOD absorption is completely anal- 
ogous to that for HOH including the two low intensity sidebands. The 
spectrum of HOD in carbon tetrachloride shown in Figure 3 for an eighteen 
to one HOD to DOD mixture is also observed to have the same basic 
structure for VY; and V- 

Since the solvents investigated are opaque in the region of the 


1. a search was made for the bending 


HOH bending vibration near 1595 cn 
frequencies of DOD and HOD. The HOD bending mode was found to be ob- 
scured in 50 mm path lengths of both carbon tetrachloride and carbon 
disulfide. For DOD in carbon disulfide it was possible to observe a 

low frequency band at 1180 cm! which was assigned as the bending vibra-~ 
tion, Y- In view of this appearance of V5 very close to the gas 


l, it appears that the bending vibration of 


phase frequency of 1178 cn” 
the water molecule is almost entirely unaffected by solution in carbon 
tetrachloride and probably therefore in other non-polar solvents. 
Hornig, using infrared data obtained from systems of dilute HOD 
in both liquid HOH and DOD, states that the shape of the HOH absorp- 
tions in the stretching regions results from overlap of Vi» V3? and 


2%; and that the intensity variations are caused by Fermi reSonance 


between 2V>5 and 1° 21,22] This type of system exhibits a decrease 


1 1 


of about 300 cm™” in ‘,, and an increase of about 50 cm™ in V> when 
compared to gas phase frequencies. These shifts bring 2 Y,, and 7; 
sufficiently near one another to allow the possibility of Fermi resonance. 
Considering the present case of HOH in carbon tetrachloride, it 
seems by analogy with DOD in carbon disulfide, that the shift of ve 
from its gas phase frequency is essentially zero. The shift of u- is 
-1 


a decrease of about 50 cm The 27, mode for HOH vapor is observed 
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at 3m on, lowered from the harmonic potential value of 3190 en. 


Considering the unshifted position of V5 in non-polar solvents it is 
reasonable to conclude that it is very unlikely that 2Y%, will exceed 
3190 em! in solution in non-polar solvents. This would eliminate 27, 
aS an explanation of the 3540 em~t sideband. These considerations lead 


to a minimum spacing between 2V> and a4 of about 400 cm! 


imp lying 
that Fermi resonance is unlikely as an explanation of the broad, low 
intensity band at 3540 em~!. This disposition of 2V5 and Fermi 
resonance leaves as the remaining plausible explanation of the band 
the possibility that it is an envelope of rotational lines. 

The persistance of rotational structure of a vibrational band of 
a molecule passing from a vapor to a liquid phase has been given little 
theoretical treatment. Longuet-Higgins expects that hindered rotation 
may persist through this transition, and that for small, nearly spheri- 
cal molecules the rotational envelope should be recognizable in favor- 
able cases. |23| The interpretation of the 3804 em-! and 3540 em7! side- 
bands as rotational envelopes is in this view a reasonable assignment. 

If this assignment is correct it might be possible to observe a 
change in these envelopes as a function of temperature. The spectrum 
of an asymmetric rotor such as water will not show the simple changes 
observed for a diatomic molecule. However, an increase in temperature 
giving an increase in the population of higher energy levels would pro- 
duce the same basic change in a gross feature such as the rotation en- 
velope; an extension of the envelope to a wider frequency range with 
greater relative intensity for the extremes of the envelope. In con- 
trast to this expected increase of intensity with temperature, the re- 


sults of Hughes et al. (that the intensity of the OH absorbance decreases 


with increasing temperature) leads to the possibility of either an in- 
crease or decrease. hs In fact, no change in intensity of the envelope 
was was identified over the 3°C to 60°C range observed. Some tentative 
changes were noted but were not reproducible. The intensities of the 
vibrational absorption maxima did decrease with increasing temperature 
in approximately the manner described by Hughes for all three species 

of water molecules in both solvents. h5| Due to the difficulty of 
accurately measuring integrated intensities, only maximum optical den- 
Sities were noted. Using the cells described, maintenance of a uniform 
temperature was not possible. It would be desirable to use an apparatus 
of the type used by Hughes: a Specially constructed metal compartment 
which maintains an accurate, uniform temperature over a wide range. hs 
Some manner of sealing the system under pressure would also be useful. 

The fact that the 3804 cm~! and 3540 em=! sidebands did not show 
a detectable decrease of intensity over this temperature range indicates 
that these bands are very probably not overtone or combination bands 
since the temperature effect would have markedly reduced their intensi- 
ties. 

It seems, then, that although there is little direct evidence that 
these two Sidebands are rotational envelopes, the dismissal of all other 
proposed sources for the bands leaves this as the most probably inter- 
pretation. 

The hygroscopic natures of carbon tetrachloride and carbon disul- 
fide are exhibited in the great difficulty encountered in drying the 
solvents and in the rapid re-absSorption of the water. A glass stoppered 
cell of carbon tetrachloride began to show HOH absorptions in periods 


on the order of one hour. Almost paradoxically, however, the solvent is 
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Saturated by a mere one part in 104 of water. ho] The tenacious grip 
of carbon tetrachloride and carbon disulfide on water may well be ex- 
plained by interstitial trapping. Such a mechanism would reconcile the 
observed difficulty in removing trace amounts of water with the low 
solubility while also allowing the water molecule to remain essentially 
Free of interactions. 

A normal coordinate analysis of the data was performed using the 
methods of Wilson, Decius, and Cross. 7 The experimental values of 
the fundamental vibration frequencies were used to calculate the quad- 
ratic force constants in the potential energy expression. The constants 
treated are those for a modified valence force treatment. In the case 
of water, there are four constants. These are for OH bond stretching, 
H-O-H bond angle deformation, and two constants for bond interactions 
(stretch-stretch and bend-stretch). That this treatment is approximate 
will readily be seen by noting that anharmonicity terms (cubic, quartic, 
etc.) are omitted from the potential energy. Calculation of precise 
quadratic force constants requires the use of "mechanical" frequencies 
of vibration, the frequencies which would be observed if there were no 
anharmonic terms in the potential energy. This calculation requires a 
substantial quantity of data regarding overtone and combination levels 
of the molecule and was accomplished for HOH vapor originally by Darling 
and Dennison. |1| 

The four valence force constants were computed here using five fre- 
quencies CV; and V3 of HOH, ae Vs and Y 3 of DOD) and alternately 
with seven frequencies, (the original five plus al and 3 for HOD). 

A best fit for each calculation was obtained using a least squares tech- 


nique. The computation was performed on a Control Data Corporation 1604 
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TABLE LII 
Results of Calculations 


All Frequencies in em! 








Species Mode Observed Calculated 
A B 
HOH yon 3614 3608 3617 
YV 3 3705 sim | 3705 
DOD VY 1 2642 2645 2646 
V5 1180 1180 1180 
Y 2 2foal ‘asa bp) Jaffa 
HOD y 2692 2704 -~ 
; 3660 3655 =~ 


Force Constants* 


Constant Bao 2 oa 
7.24 4.67 8.43 

0.92 Seeoe Oo 77 

0.07 -2.50 -0.10 

-0.68 ~0.84 0.25 


*Units: Millidynes per Angstrom 


Calculation A: Seven frequencies input (HOH, DOD, and HOD) 
Calculation B: Five frequencies input (HOH and DOD only) 


Calculation C: Given by Wilson, Decius and Cross as the mechanical 
frequencies for HOH. 7) 
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digital computer using FORTRAN 60 input. Details of the calculation are 
presented in the appendix. Results are presented in Table III. 

All of the factors discussed above indicate the conclusion that the 
water molecule when dissolved in a non-polar solvent Such as carbon tet- 
rachloride or carbon disulfide interacts very weakly with the solvent. 
In summary these factors are: 

(a) It has been shown that it iS consiStent with all observations 


: and 3540 em~! bands as envelopes of rotation 


to regard the 3804 cm 
lines, indicating that the molecules are relatively free to rotate; 

(b) The very low frequency shifts observed upon change from vapor 
to condensed phase indicates that the potential field of the molecule 
is nearly undisturbed by molecular association; 

(c) A calculation based on an unaSsociated molecule is reasonably 
consistent; 

(d) The hygroscopicity of carbon tetrachloride and carbon disul- 
fide is related to the low Saturation value of water by a mechanism of 


interstitial trapping which is consistent with minimum molecular inter- 


actions. 
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APPENDIX 


The quadratic force constants for the potential energy expression 
for the water molecule were calculated by fitting the calculated vibra- 
tional frequencies to the observed frequencies by the method of least 
Squares. Calculated vibrational frequencies were obtained through 
solution of the Secular equation in the form given by Wilson, Decius 
and Cross, 

GE-E| =0 


where G, E, and F are matrices. h7| G is the Wilson G matrix given by 


H, +43 M4, cos G —— 
G =] 43 cos G My + My -M, sin G 
-H 2 sin @ -4» sin @ - a 
— = 7: [4 1 +y*2M,(1-cos 2) | 


for non-linear triatomic molecules such as HOH, DOD and HOD. F is 
the matrix of quadratic force constants. For water these constants 
are f, bond stretching; fx ; bond angle bending; fer stretching inter- 


action; and E nok » Sstretch-bend interaction. 


aa trary a 
ayy = tr ory Ero * Tk 
f f e4 


Tic TOK 


E is the third order identity matrix. The machine solution of the 
equations is quite simple: enter with values for G and F, perform 
the indicated product and then diagonalize this product. The resulting 


eigenvalues are the desired fundamental vibration frequencies. 


“ry — 


The least square fit of data was derived as follows: 


Solve the system of equations: 


3 a bo 
Sap Reuss" 


where g (25) are calculated frequencies, f; are observed frequencies 


and Z are the force constants. Differentiating: 


a 8} 


= y 2(f£. - g-) —=———— = () 
yi 
i ows * 
Expanding g. in terms of 233 
= eo > &i : : 
S > 8. ia —— A Zy neglecting higher terms so that 


k S 2k 


g; oe 
Mi [i -e- Lh Sy-44] Sao 


1 

Let 

9 8i O 

S72, eee Zp se and €, #f, - gj 

Then 

2 Ee, 2 y, iA | Ds = 6 

1 k 

cone 2h Ds, 4, D.. = 0 

Lk 

7 i ee ed 

d, me eae >, D.. D..A, = 90 

i jl L ik g5 Lk k 
Thus 


: CD BiAonee= 0 
(Dé); ( 4d), 


where now 
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(D € ds and (D DA) 


are column matrices. 


(D'D) is a square matrix, hopefully with an inverse. Then 


I 


(Di D) ae Seno 


So that 
A =a De 


For the second iteration new force constants (25) are 


awe = Zanes where all are column matrices. 
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